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^ ■ Abstract 
O 

| 7j i We calculate the nuclear modification factor for RHIC and LHC conditions accounting 
ly-^ ■ for the radiative and collisional parton energy loss with the running coupling constant. 

We find that the RHIC data can be explained both in the scenario with the chemically 
^ . equilibrium quark-gluon plasma and purely gluonic plasma with slightly different thermal 
^! suppression of the coupling constant. The role of the parton energy gain due to gluon 
Q-i| absorption is also investigated. Our results show that the energy gain gives negligible 
■ effect. 

(N 

^ ! 1. It is widely believed that suppression of the high-pr hadrons in AA-collisions (jet 

I quenching (JQ)) observed at RHIC (for a review, see [1]) is dominated by the induced 
^ I gluon emission [2, 3, 4, 5, 6, 7] in the hot quark-gluon plasma (QGP) produced at the 
initial stage of AA-collisions. There are currently considerable theoretical efforts in the 
development of the quantitative methods for computation of JQ [8, 9, 10, 11, 12, 13] 
OO ! which can be used for the tomographic analysis of the QGP. In the present paper we 
study JQ using the light-cone path integral (LCPI) approach to the radiative energy loss 
[3, 4]. In this formalism the probability of gluon emission is expressed through the Green's 
^ ■ function of a two-dimensional Schrodinger equation with an imaginary potential. This 
approach has not the restrictions on the applicability of the BDMPS approach [2] (valid 
only for massless partons in the limit of strong Landau-Pomeranchuk-Migdal effect) and 
the GLV formalism [6] (applicable only to a thin plasma in the regime of small Landau- 
Pomeranchuk-Migdal suppression). We perform the calculations with accurate treatment 
of the Coulomb effects. If one neglects these effects the gluon spectrum can be expressed 
in terms of the oscillator Green's function and the medium may be characterized by the 
well-known transport coefficient q [2, 4]. However, the oscillator approximation can lead to 
uncontrolled errors since it gives a physically absurd prediction that for massless partons 
the dominating = 1 rescattering contribution vanishes [14, 15]. Besides the radiation 
energy loss we include the collisional energy loss. Both the contributions are calculated 
with the running coupling constant. Also, we investigate the impact of the parton energy 
gain due to gluon absorption from the QGP on JQ. 

We calculate the nuclear modification factor Raa, which characterizes JQ, accounting 
for the fluctuations of the parton path lengths in the QGP. In the treatment of multiple 
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gluon emission we use a new method which takes into account time ordering of the DGLAP 
and the induced radiation stages. We compare the theoretical results with the data 
obtained at RHIC by the PHENIX Collaboration[16] and give prediction for LHC. Our 
principle purpose in comparing with the RHIC data is to understand whether the observed 
JQ is consistent with the entropy of the QGP required by the hydrodynamical simulations 
of the AA-coUisions for reproducing the observed particle multiphcities. Our results show 
that JQ and particle multiplicities can be naturally reconciled. Contrary to the conclusion 
of Ref. [17] that the observed at RHIC JQ is consistent only with purely gluonic plasma, 
we find that the scenario with the chemically equilibrium QGP is also possible. A good 
description of the JQ RHIC data can be obtained in this scenario with the thermal 
suppression of the coupling constant qualitatively consistent with the lattice results. 

2. As usual we define the nuclear modification factor for 7474-coUisions as 

^ dN{A + A^h + X)/dpTdy 
^ TAA{b)da{N + N^h + X)/dpTdy ' ^ ^ 

where p^- is the hadron transverse momentum, y is rapidity (we consider the central 
region y = 0), b is the impact parameter, TAA{i>) = J dpTA{p)TA{p — h) , Ta is the nucleus 
profile function. The differential yield for high-p^ hadron production in A^l-coUision can 
be written in the form 

dprdy J dprdy 

where dam{N + N ^ h + X)/dpTdy is the medium-modified cross section for the N + 
N ^ h + X process. In analogy to the ordinary pQCD formula, we write it in the form 

damjN + N ^h + X) _^ f^d^r^m r da{N + N^i + X) 
d^ -^Jo ^^V^^^'^^ d^ ■ 

Here = pr/^ is the parton transverse momentum, D^^- is the medium-modified frag- 
mentation function (FF) for transition of the parton i to the observed hadron h, and 
da(N + N ^ i + X)/dplpdy is the ordinary hard cross section. For the parton virtuality 
scale Q we take the parton transverse momentum p^. We assume that hadronization of 
the fast partons occurs after escaping from the QGP. This hadronization process should 
be described by the FFs at relatively small fragmentation scale, fih- Indeed, from the un- 
certainty relation AE'At ^ 1 one can obtain for the L dependence of the parton virtuality 
Q^{L) ~ max (Q/L, (5o), where we have introduced some minimal nonperturbative scale 
Qo ~ 1 ~ 2 GeV. For RHIC and LHC conditions the size of the QGP is quite large (> Ra, 
where Ra is the nucleus radius), and from the above formula one sees that for partons 
with energy E < 100 GeV the hadronization of the final partons may be described by the 
FFs at the scale ^ Qo- Then we can write 

/■I dz' 

D^llz, Q) ^ —Dn,Az/z\ Qo)D]},{z', Qo, Q) , (4) 

where Dh/j{z,QQ) is the FF in vacuum, and D^^-^z' ,Qq,Q) is the medium-modified FF 
for transition of the initial parton i with virtuality Q to the parton j with the virtuality 
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Qo- Presently there is no a systematic method for calculation of the medium- modified 
FFs which treats on an even footing the DGLAP and induced radiation processes. In 
the present paper we use the picture based on the time ordering of the DGLAP and 
the induced radiation stages which should be a reasonable approximation for not very 
high parton energies, say, E < 100 GeV. It uses the fact that at such energies the typical 
length/time scale of the DGLAP stage is smaller than the longitudinal scale of the induced 
radiation stage. The gluon emission scale for the DGLAP stage can be estimated using the 
gluon formation length lp{x, k^) ~ 2Ex{l — x)/ (k^ + e^), where x is the gluon fractional 
longitudinal momentum, and e in terms of the effective parton masses reads = rrigx'^ + 
m^{l — x). Using the vacuum spectrum of the gluon emission from a quark 



dN Cpasik^) f 2/o^^ ^'^ 



1- + -V2))vr^ (5) 



dkj^dx nx ^ ' {k\ + e 

one can obtain for the typical formation length Ip ~ 0.3 — 1 fm for E < 100 GeV (if 
one takes ruq ~ 0.3 GeV and rrig ~ 0.75 GeV [18]). This estimate is obtained in the 
one gluon approximation. However, it should be qualitatively correct since in the energy 
interval of interest the number emitted gluons is small Ng <2, and the first hardest gluon 
dominates the DGLAP energy loss. Thus we see that the DGLAP time scale is about the 
formation time for the QGP, tq ~ 0.5 — 1 fm. Since the induced radiation is dominated 
by the distances from L ~ Tq up to L ~ Rj^ one can neglect the interference between the 
DGLAP and the induced radiation stages. In this approximation we can write 

DT/iiz, Qo, Q) = I ^D'l/fiz/z', E{)Df,f^^^{z\ Qo, Q) , (6) 

where Ei = Qz', DJjf is the induced radiation FF (it depends on the parton energy E, but 
not the virtuality), and Dfjf^^"^^ is the DGLAP partonic FF. In numerical calculations 
the DGLAP FFs have been evaluated with the help of the PYTHIA event generator [19]. 

The induced radiation FFs have been calculated making use the probability distribu- 
tion of the 1 — s> 2 partonic processes obtained in the LCPI approach. We have taken into 
account only the processes with gluon emission, and the process g ^ qq which gives a 
small contribution has been neglected. For calculation the one gluon emission distribution 
we use the method elaborated in [20] . To calculate the DJjf one needs to take into account 
the multiple gluon emission. Unfortunately, up to now, there is no an accurate method 
of incorporating the multiple gluon emission. We follow the analysis [8] and employ the 
Landau method developed originally for the soft photon emission. In this approximation 
the quark energy loss distribution has the form 
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dPi 



duj 



6 AE — ^-^P 



=1 



, dP 
dui—— 
duj 



(7) 



where dP/duj is the probability distribution for one gluon emission. This approxima- 
tion leads to the leakage of the probability to the unphysical region of AE > E [9]. 
To avoid the quark charge non- conservation we define the renormalized distribution 
P{AE) = KgP{AE) with Kg ^ dAEP{AE)/J^ dAEP{AE). We use the renor- 
malized distribution to define the in-medium FF L)*"^(2;) = P(A£' = £'(1 — z)). To ensure 
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the momentum conservation we take into account the q ^ g transition as well. At the one 
gluon level the corresponding FF can be written as D^J^^'^{z) — dP{uj = zE)/duj. This auto- 
matically leads to the FFs which satisfy the momentum sum rule. We use the same form 
of the q ^ g FF for the case with the multiple gluon emission. To satisfy the momentum 
sum rule (which are not valid after the renormalization of the q ^ q distribution) we 
multiply it by a renormalization coefficient Kg defined from the total momentum conser- 
vation. This procedure seems to be reasonable since the nuclear modification factor are 
only sensitive the behavior of the FFs at z close to unity [8] where the form of the q ^ g 
distribution should not be very sensitive to the multiple gluon emission. In the case of 
the g ^ g transition we use the following prescription. In the first step we define 1^*"^ at 
z > 0.5 through the Landau distribution P{/S.E), and in the soft region z < 0.5 (where the 
multiple gluon emission and the Sudakov suppression strongly compensate each other) we 
use the one gluon distribution. Then we multiply this FF by a renormalization coefficient 
Kg to ensure the momentum conservation (since the number of gluons is not conserved 
the arguments based on the conservation of the probability cannot be used in this case) . 

In the above discussion we ignored the coUisional energy loss. Presently there is 
no an accurate method for incorporating of the coUisional energy loss in the scheme of 
the medium-modified FFs. In the present work we view the coUisional energy loss as a 
perturbation and incorporate it into our model by a small renormalization of the QGP 
density according to the change in the A£' due to the coUisional energy loss. To evaluate 
the coUisional energy loss we use the Bjorken method [21] with an accurate treatment 
of kinematics of the binary collisions (the details can be found in [22]). We use the 
same infrared cutoffs and parametrization of the coupling constant for the radiative and 
coUisional energy loss, which is important for minimizing the theoretical uncertainties in 
the fraction of the coUisional contribution. 

We calculate the cross sections for the N + N ^ q{g) + X processes using the LO 
pQCD formula with the CTEQ6 [23] parton distribution functions. To account for the 
nuclear modification of the parton densities (which leads to some small deviation of Raa 
from unity even without parton energy loss) we include the EKS98 correction [24]. To 
simulate the higher order K-factor we follow the prescription used in the PYTHIA event 
generator [19] with replacement of the Q argument of ag by a lower value cQ. We take 
c = 0.265 which allows to describe well the data on tt*^ p^-spectrum in the pp-coUisions. 
For the FFs Dh/q(g){z, Qq) we use the KKP parametrization [25]. 

3. For calculations of the induced gluon spectrum and coUisional energy loss with the help 
of the formulas given in [20, 22] we must specify the form of the coupling constant and 
the mass parameters (the quasiparticle masses and the Debye mass). In our calculations 
we use the running coupling constant. We parametrize as{Q^) by the one-loop expression 
and assume that it is frozen at some value q;{^ for Q < Qf^.. This form with 0.7 
(Q/r a; 0.82 GcV for Aqcd = 0.3 GeV) allows one to describe well the HERA data on 
the low- a; structure functions within the dipole approach [26, 18, 27]. A similar value 

of af follows from the relation ^"^^dQ^^^ ^ 0.36 GeV obtained in [28] from the 
analysis of the heavy quark energy loss in vacuum. In vacuum the stopping of the growth 
of as at low Q may be caused by the nonperturbative effects [28]. In the QGP thermal 
partons can give additional suppression of as at low momenta (Q ~ 2 — ST). The lattice 
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simulations [29] give as{T) smoothly decreasing from ~ 0.5 at T 175 MeV to ~ 0.35 
at T a; 400 MeV. However, the thermal as{T) in some sense gives the mean value of a^. 
For this reason one can expect that the thermal as{T) should be somewhat smaller than 
the in-medium a{^. To clear up whether the RHIC data on jet quenching agree with the 
thermal suppression of ag we perform numerical calculations for different values of 

As in [22] we use the quasiparticle masses obtained in Ref. [30] from the analysis 
of the lattice data within the quasiparticle model. For the relevant range of the plasma 
temperature T ~ (1 — 3)Tc the analysis [30] g IVGS 777-/7 0.3 and rUg ^ 0.4 GeV. To fix 
the Debye mass in the QGP we use the results of the lattice calculations for Nf = 2 [31] 
which give the ratio Hd/T slowly decreasing with T (^d/T ^ 3 at T 1.5Tc, /xd/T 2.4 
at T~ 4Te). 

4. We describe the QGP in the Bjorken model [32] with the longitudinal expansion 
which gives the proper time dependence of the plasma temperature T^r = TqTq (Tq is 
the initial plasma temperature). To simplify the numerical calculations for each value 
of the impact parameter b we neglect the variation of Tq in the transverse directions. 
For each b we define its own effective initial temperature evaluated with the help of 
the entropy distribution adjusted in the hydrodynamic analysis [33] of the RHIC data 
on the small-px hadron spectra. For Au+Au collisions at ^/i = 200 GeV it reads [33] 
dS{T, p, h)/dpdz = ^(^^^^ [cidNpart/ dp + (1 — a)dNcoii/ dp] , where C = 24, and a = 0.85, 
Npart and Ncou are the number of participants and binary collisions evaluated in the 
Glauber model. In evaluating the entropy distribution we use the Woods-Saxon nucleus 
density pA{r) =Cnorm/ + exp[(r — c)/d]} with c = 1.07A^/^ fm, and d = 0.545 fm. For 
central Au+Au collisions at ^/s = 200 GeV this gives To ~ 320 MeV for tq = 0.5 fm. 
It was assumed that the QGP occupies the region r < (c + kd), with k = 2 (for k = 1 
Raa changes slightly). To perform the extrapolation to the LHC energy y/s — 5500 GeV 
we use the energy dependence of the entropy similar to the energy dependence of the 
total particle rapidity density oc iVpart In (v^/1.5) observed at RHIC [34]. It gives for 
central Pb+Pb collisions To ~ 404 MeV. The fast parton path length in the QGP, L, 
in the medium has been calculated according to the position of the hard reaction in the 
impact parameter plane. To take into account the fact at times about 1 — 2 units of Ra 
the transverse expansion should lead to fast cooling of the hot QCD matter [32] we also 
impose the condition L < L^ax- We performed the calculations for two values L^ax — 6 
and 8 fm. 

5. We present the numerical results for a('^ = 0.7, 0.5 and 0.4. The higher value seems 
to be reasonable in the absence of the thermal effects since it comes from the analyses 
of the HERA data on the low-x structure functions [18] and heavy quark energy loss in 
vacuum [28]. In Fig. 1 we plot Raa for 7r° production in the central Au-|-Au coUisions at 
^/s — 200 GeV. The theoretical curves corresponds to L^ax — 8 fm. The choice L^ax — 6 
fm gives Raa higher just by about 3-8%. The experimental points in Fig. 1 are from [16]. 
The upper and lower panels show the results for the chemically equilibrium and purely 
gluonic plasmas, respectively. The results are presented for the purely radiative energy 
loss and with inclusion of the coUisional energy loss and the radiative energy gain. We 
have found that the effect of the radiative energy gain on Raa is practically negligible 
and can be safely neglected. Besides the total (quarks plus gluons) contribution in Fig. 1 
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we also show separately the contributions from quarks and gluons. The grows of Raa for 
gluons is due to the q ^ g transition which is usually neglected. However, it does not 
affect strongly the total Raa since for ^/s = 200 GeV the gluon contribution to the hard 
cross section is small at px ^ 15 GeV. As can be seen from Fig. 1 the coUisional energy 
loss suppresses Raa only by about 15-25%. This is in contradiction with strong coUisional 
JQ found in [35] . To illustrate the effect of the time ordering of the DGLAP and induced 
radiation stages in Fig. 1 we also plot the results for the inverse order of these stages. 
One can see that the results for these two prescription are very close. This fact may be 
explained by the dominance of the soft gluon emission at RHIC energies. 

Prom Fig. 1 one can see that the theoretical Raa for the chemically equihbrium plasma 
obtained with q;{^ = 0.5 is in qualitative agreement with the experimental one. The 
scenario with purely gluonic plasma can be consistent with the RHIC data if ~ 0.4. 
In the light of the lattice results [29] the value ~ 0.5 seems to be reasonable for the 
RHIC (and LHC) conditions. Thus one sees that, contrary to the conclusion of Ref. [17], 
the scenario with the chemically equilibrated plasma can not be excluded. Note that the 
distribution of the entropy [33] used in our calculations gives the entropy rapidity density 
dS/dy ~ 6500. If we take a smaller value dS/dy ^ 5100 obtained in [36] approximately 
the same Raa as for a^^ — 0.4 and 0.5 can be obtained with the values a^^ ~ 0.45 and 
0.55. 

In Fig. 2 we show Raa as a function of Npart- One sees that the model reproduces 
qualitatively the growth of Raa with decrease of Npart- But for very peripheral collisions 
with small Npart it overestimates the observed Raa- This may be connected with inade- 
quacy of the neglect of the transverse motion of the matter for thin plasma. Also, in this 
region the neglect of the variation of Tq in the impact parameter space may be inadequate 
as well. Probably for similar reasons the model underestimate the ellipticity parameter V2 
(which is also sensitive to the evolution of the QGP for the peripheral coUisions) . Our cal- 
culations give V2 ~ 0.05 — 0.08 for ~ 5 GeV while the experiment gives V2 ~ 0.1 — 0.15 
[37]. 

In Fig. 3 we plot the theoretical results similar to that shown in Fig. 1 but for Pb+Pb 
collisions at LHC for ^/s = 5500 GeV. One can see that the effect of the coUisional 
energy loss becomes smaller for LHC conditions. But the effect of the time ordering 
of the DGLAP and induced radiation stages is bigger as compared to the RHIC. The 
difference between the results for the chemically equilibrium and non-equilibrium plasmas 
is relatively small. 

6. In summary, we have calculated the nuclear modification factor for the RHIC and 
LHC conditions accounting for both the radiative and coUisional energy losses with the 
running ag- The radiative energy loss has been calculated within the LCPI approach 
[3] . The coUisional energy loss has been evaluated in the Bjorken model of elastic binary 
collisions with an accurate treatment of kinematics of the binary collisions. In contrast to 
[35] we find relatively small effect of the coUisional energy loss on JQ. We also investigated 
the effect of the parton energy gain due to the induced gluon absorption. We find that 
this effect is negligible for the RHIC and LHC conditions. 

The calculations are performed using a new algorithm for the multiple gluon emission 
which takes into account the time ordering of the DGLAP and induced gluon emission 
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stages. We find that the effect of the ordering of these two stages is relatively small for 
RHIC conditions, but becomes bigger for LHC. 

Comparison of our theoretical results with the RHIC data show that Raa can be 
described in the scenario with the chemically equilibrium QGP with the entropy extracted 
from the hydrodynamical simulation of the AA coUisions at RHIC energies. The scenario 
with the purely gluonic plasma is also possible, but requires somewhat stronger thermal 
suppression of a^. This contradicts to the conclusion of Ref. [17] that the observed JQ 
and total entropy of the QGP are incompatible with the chemically equilibrium plasma 
scenario. 
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Figure 1: The nuclear modification factor Raa for tt" production in the central Au+Au 
collisions at = 200 GeV for a^/ = 0.7 ((a),(d)), = 0.5 ((b),(e)), a^/ = 0.4 ((c),(f)). 
The upper panels are for the chemically equilibrium plasma, and the lower ones for purely 
gluonic plasma with the same entropy. For all the theoretical curves L^ax = 8 fm. Solid 
line: the total (quarks plus gluons) radiative Raa- Dashed line: the radiative Raa for 7^° 
from quarks. Long-dash fine: the radiative Raa for ti"" from gluons. Dash-dotted line: the 
total (quarks plus gluons) Raa including the radiative plus coUisional energy loss and the 
energy gain due to gluon absorption. Dotted curves show the total (quarks plus gluons) 
radiative Raa for the inverse time order of the DGLAP and induced gluon emission stages. 
The experimental points are the data obtained by the PHENIX Collaboration [16] for the 
most central (0-5%) coUisions. 
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Figure 2: The nuclear modification factor Raa for ti^ production in Au+Au at ^/s = 200 
GeV for Pt > 5 GeV (left panels) and pt > 10 GeV (right panels) as a function of Npart- 
The upper panels are for the chemically equilibrium plasma, and the lower ones for purely 
gluonic plasma with the same entropy. The theoretical curves show the total (quarks plus 
gluons) Raa including the radiative plus coUisional energy loss and the energy gain due 
to gluon absorption for = 0.7 (solid line), a{'^ — 0.5 (dashed hne), a{'^ — 0.7 (dotted 
line). The experimental points are from [16]. 
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